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ABSTRACT 



The hydraulic analogy to the isentropic flow of a 
compressible fluid has long been recognized.. This study was 
undertaken to find some means, if possible, of extending the 
analogy so as to simulate the addition of heat in compressible 
fluid flow. It was shown that while the analogy does not 
hold for this phenomenon or for any non-isentropic process, 
a method does exist by which the isentropic addition of energy 
may be simulated. 

A secondary aim of the investigation was to study and 
improve the performance of the laboratory water table of the 
Department of Aeronautics, United States Naval Postgraduate 
School. The objective was to make the equipment suitable 

T 

both for flow visualization demonstrations and for further 
research purposes. 

The water table was rebuilt to include a lighted surface 
incorporating a two-inch grid for photographic reference. An 
accurate method of measuring local and reservoir water depths 
was established. 

Photographs were taken of flow about models of various 
shapes to demonstrate the usefulness of the water table and 
the hydraulic analogy in the study of such flows. 
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1. Introduction 



The hydraulic analogy has long been recognized as a 
qualitative and, to a limited extent, quantitative aid in 
the study of compressible fluid flow. [2] Its basis is 
the similarity of the equations governing an incompressible 
fluid flowing in a channel with a free surface and those 
which describe the flow of a compressible gas. Its useful- 
ness lies in the fact that supersonic flow may be simulated 
about various shapes for prolqnged periods of time at very 
low cost. It also allows the operator to make numerous var- 
iations in the flow patterns quite rapidly and inexpensively. 

It was desired to find a method by which the analogy 
could be extended to include the addition of energy to the 
stream. Results analogous to those obtained by heating a gas 
flowing in a constant area duct (Rayleigh line) were expected 
and, it was thought, a phenomenon analogous to choked flow 
might develop. 

Coincident with the above investigation, the water table 
of the Department of Aeronautics, United States Naval Post- 
graduate School, was to be refurbished for use as a teaching 
aid and as a tool for later investigations of various super- 
sonic flow phenomena. 

2. The Hydraulic Analogy 

The hydraulic analogy to two-dimensional frictionless 
isentropic fluid flow has its basis in the similarity of the 
equations of motion governing the compressible flow situation 
and those governing the quasi-two-dimensional flow of an in- 
compressible fluid in a horizontal channel with a free surface. 
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The latter is called quasi- two-dimensional since there are 
depth changes from point to point, but vertical fluid velo- 
cities and accelerations are small and are neglected. [1] 
The energy equations governing the two flows are ex- 
amined below. 



Compressible Fluid 
h + V 2 /^ + gz = h Q = const. 
Assuming horizontal flow, 
z = a constant. Solving 



for V, 



V =V2(h Q -h) 



Incompressible Fluid 
p/pg + V 2 /2g + y = y 0 = const. 
Assuming p constant at all 
points and solving for V, 



V = V2g(y 0 -y) 



The continuity equations for the two flows are also 
analogous. This has been shown to be true for two dimensional 
flow [ 4 "| , but for simplicity it is shown below for the one- 

dimensional case. 

Incompressible Fluid 
pVw y = m = constant 



Compressible Fluid 
pAV = m = constant 



Letting h = CpT, and re-writing the energy equation for 
a compressible fluid in terms of T, k, and M, it is seen that 



T 0 /T = 1 + (M) 2 

The result of re-writing the energy equation for the 
compressible fluid is, in terms of y and F, 

y 0 /y = 1 + £(F) 2 

It may now be seen that the analogy holds If the compress- 
ible fluid is such that k = 2.0. This value of k is peculiar 
to a channel of rectangular cross section. Alternative cross 
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sections allow other values of k. A triangular cross sec= 
tion, for instances, corresponds to k = 1«5» Use of these 
alternative cross sections has been discussed by Shapiro [71 
and [9] . 

2 

Pressure ratio, p/p q9 is analogued by (y/y 0 ) . This 
is shown by an examination of the equation of state, p = pRT. 

p 

Since y is the analogue of both p and T, y must be analogous 
to p. 

The expression for wave propagation speed (sonic velocity) 
in a compressible fluid iss 

c = VKRT 

The wave propagation speed for the incompressible sit” 
uation iss 

a = Vgy 

This expression is not valid in the ocean or in deep 
channels, since the assumption of large wave lengths compared 
to depth is explicit in its deviation. [41 

Thus dimensionless velocity rations may be written for 
both cases. 

Compressible Fluid Incompressible Fluid 

M = V/c F = V/a 

It is of Interest to note that abrupt transition from 
flow having a Fronde number greater than one (supercritical 
or shooting flow) to a flow having a Froude number less than 
one ( subcritical, tranquil, or streaming flow) is accompanied 
by a hydraulic jump. This is somewhat analogous to the shock 
wave which accompanies sudden transition from supersonic to 
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subsonic flow. It is not completely analogous, since 
the stagnation temperature of a compressible gas remains 
constant through the shock, while the total head y Q de- 
creases across a jump due to the dissipation of mechani- 
cal energy in the turbulence of the jump. Figure 1 com- 
pares the conditions prevailing on both sides of a shock 
wave with those on both sides of the hydraulic jump. 



Normal Shock 


Hydraulic Jump 
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FIGURE 1 

ANALOGY BETWEEN NORMAL SHOCK AND HYDRAULIC JUMP 

Although the momentum. equation has not been used in 
the establishment of the analogy, it may be employed 
along with the equations of energy and continuity to 
determine the decrease in y Q across the jump. Assuming 
the hydrostatic pressure to vary linearly with depth the 
average pressure on a face perpendicular to the direction 
of flow will be pgy/2. The momentum equation can now be 
written (assuming negligible external friction): 



A 



PSy yfy p 

■ — ^ + PwyV^ = constant 

Expressing velocity in terms of the Proude number, a 2 F 2 , 

where a 2 = gy. Making these substitutions and combining the 
constant terms p, g 9 and w, the equation can be re~wrltteng 

2 2 

y (1 + 2F ) = constant 

The continuity equation can also be expressed in terms of 
y and P; 

yip = constant 

Raising the continuity equation to the second power and 
dividing by the momentum equation to the power yeilds the 



expressions ^2 

(1+2F 2 ) 

This expression is double valued except when P 



. 2 / x - constant 



1.0 The e= 

quation must hold across the hydraulic Jump, since the only 
restriction imposed upon it is that wall friction be neglig* 
lble. The equation can best be solved graphically by plot” 
ting the function P^/(l+P^)2 against F. Two values thus ob= 



tained, for example , are P = 3-0 and 0.41. See figure 2, 



The energy equation may now be used to find y 



0 / 



Since 



total energy is not constant across the jumpy the energy 
equation must be written twice. 

n . .. 



= y l + 2g 



. 4 

y 2 + 2 s 



Substituting = gyP^ and dividing, it is seen that 



y o 2 _ y 2 (i+i|/2)^ 

S = (l+sf/2? ' 

Since, from the continuity equation, yp/V^ = (P^/F^) 2 ^, the 
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ratio of the energy equations can be written: 



P. 



2/3 



l+Po/2 



y o 1 *2 \ 1 + 1^12 

Substituting the values of and Fg previously deter- 
mined (3-0 and 0.4l) y„ was calculated to be 0.74 y_ , 

°2 °1 

a decrease of 26 %, 

It is seen that the decrease in y Q across an hy- 
draulic jump is considerable at high Froude numbers, and 
of interest to note that a method of establishing this 
decrease was not found in the literature. The momentum 
equation was not applied _:or even mentioned, although it 
is satisfied and is analogous to the momentum equation 
of compressible fluid flow. 




COMBINATION OF MOMENTUM AND CONTINUITY 
EQUATIONS ACROSS AN HYDRAULIC JUMP 
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The literature refers to the "hydraulic analogy to 
isentropic compressible fluid flow", but does not offer a 
proof that this restriction must be imposed. 

Aa. examination of the laws of thermodynamics proves that 
the analogy is, restricted to isentropic processes. The com- 
bined first and second laws can be written 

Tds = dh- vdp 

Utilizing the relationship dh = c^dT, the combined first and 
second laws can be re-written: 

ds = c p - vdp/T 

Substituting v/t = R/p from the equation of state and divid- 
ing through by R 

ds c p dT _ dp 

R “ R T ~ p 

Since R = c p -cv, and c p/ cv = c /R = k^l * Makin 6 these 
substitutions, 

ds _ k dT dp 

R k-1 T p 

The quantities and ^ are analogued by and 2y ^ 

y 7 

respectively. 

d. s 

Substituting into the expression for jp , 

d§. - 2 d z - 2 d Z z o 
R " 2 y 2 y “ °* 

Thus it is seen that the analogy holds only for isentropic 
processes. 
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Table I summarizes the analogous quantities 



TABLE I 

ANALOGOUS QUANTITIES BETWEEN A 
COMPRESSIBLE GAS AND AN INCOMPRESSIBLE FLUID 



Compressible gas 


Incompressible Fluid 


i/i 0 


y/y 0 


>/p 0 


y/y 0 


p/p 0 


(y/'y 0 ) 2 


V = V2c p (T Q “T) 


v = V2g(y 0 -yj 


c = "VkRT 


a = Vg/ 


m = PAV 


A = p wy V 


M = v/c 


F = v/a 


shock wave 


hydraulic .lump 



3. Addition of Energy- 

Energy is usually added to a compressible gas flawing 
in a duct by burning fuel in the stream 9 thus increasing the 
total temperature. This process is not isentropic and cannot 
be simulated via the hydraulic analogy. The results of an 
attempt to "force" an analogy are discribed below. 

Energy can be added to the water at any point in the 
stream by dropping the level of the surface upon which it 
is flowing. This causes an effective increase in y Q with a 
corresponding change in velocity and Froude number. 

The energy and continuity equations can be combined to 
give an equation relating local depth y and stagnation depth 
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y Q at any fixed flow rate and channel width. 
The energy equation is written 

V 2 = 2g(y 0 -y) 

and the continuity equation 





pwyV = m 


or 


2. r 2 / 1^2 
y V = (pw) 


Substituting 


p 

for v from the energy equation 

y 2 [2g(y 0 -y) ] = ("£?) 


Let 


_ 3 JL (Jk ) 2 

5 “^g ' pw 


then 


2 3 3 

. y y 0 ~y = § 



This relationship is shown in Figure 3. 

The two branches of the curve shown in Figure 3 rep 
resent flow at subcritical and supercritical Froude num- 
bers. Flow occuring at values of y and y Q corresponding 
to the "nose" of the curve is at a Froude number of 1.0. 
Under these conditions y = 2/3 y • 
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FIGURE 3 



. COMBINATION OF ENERGY AND CONTINUITY 
EQUATIONS FOR AN INCOMPRESSIBLE FLUID 
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1 — — I 1 1 1 

l.o 2.0 3.0 4.0 

F 

FIGURE 4 s 

y/y o VERSUS FROUDE NUMBER 

Figure 4 is a plot of y/y Q against Froude number, 
based on the non-dimensialized energy .relationship, 

l+F 2 /2 

The supercritical branch of Figure 3 shows there 
is no analogy to the "choked" condition which occurs 
when heat is added to a compressible gas flowing in a 
constant area duct. As y Q increases, y decreases and 
Froude number increases without limit. In the non-' 
isentropic addition of heat to a compressible gas, T 
increases as T Q is increased and Mach number approaches 
1.0 as a limit. Choking occurs when the product pY is 
a maximum, at which point M 1.0. Since the analogue of 




pV is Vy which remains constants there will be no analogue 

to choked flow. 

The subcritical branch indicates that as y is increase 
ed y also increases. This can be visualised by thinking of 
water in a shallow stream running over a slight drop with 
little or no apparent change in water level. Here y is in- 
creased and so is y. Velocity decreases as dees Fronde number. 

If the drop is made great enough,, however, there will be 
a decrease in y, and the velocity after the drop becomes great” 
er than the velocity before it. The Froude number also in” 
creases from subcritical to supercritical. This must cor- 
respond to a shift from point A to point B as shown in Figure 
3 since flow at these points can occur en either branch at the 
same y » The magnitude of the increase in y required to pro- 
duce this shift is thought to be a function of downstream 
conditions such as channel friction and is a subject for 
further investigation. It was not carried out at this time 
since to do so would have destroyed the table surface. 

Another method of investigating the effect of an increase 
in y Q on the subcritical branch was attempted. It consisted 
of forcing a.n increase in y, by lowering the sluice and mea~ 
suring the resultant increase in y. 

No increase in y was noted as the sluice was lowered, 
even though the initial sluice position was clear of the flow. 
This was attempted at flow rates varying from 15 to 30 gallons 
per minute. These flow rates caused Fronde numbers from 0.4 to 
0.65 with the table level and the sluice in the clear position. 

It Is felt that the effect expected was not observed due 

to several factors, the foremost being the limited capacity of 
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the pump. Operating at a flow rate of 30 gallons per minute 
with the sluice lifted clear of the flow, a velocity of 0.68 
feet per second was achieved with the table level. This cor- 
responds to a reservoir depth (y ) of 0.488 inches and a 
stream depth (y) of 0.406 inches. Even though the difference 
of 0.088 inches is measurable with the apparatus described 
in Section 4, any increase in y caused by contact of the sluice 
with the water surface was not. Any further depression of the 
sluice caused a shift to supercritical flow immediately down- 
stream of the sluice. A hydraulic Jump then occured immed- 
iately and the flow was back to subcritic al. 

It is felt that if larger flow rates were possible the 

increase in the quantity (y Q -y) for subcritical flows would 

allow the demonstration of the increase in y with y described 

o 

by the equations of motion and continuity. 

This increase is shown diagramatically in Figure 5 as 
it is expected to appear. The flow in Figure 5 (a) is sub- 
critical and the sluice is set at d-^. As the sluice is de- 
pressed, in Figure 5 (h), d decreases to dg and y increases 
to yg» The Froude number also decreases. As the sluice is 
depressed further to some setting d^> transition to super- 
critical flow occurs at y 0 ^ and y^ with a Froude number 
greater than 1.0 as shown in Figures 5 ( c ) and 5 (d). 

Another method of increasing the energy of the stream 
would be to add water at some point in the stream. If this 
water were at a higher head than the original y Q the resultant 
flow would have a higher head per unit mass than the original 
stream. 
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FIGURE 5 

\ / 

FLOW AT VARIOUS SLUICE SETTINGS 

To determine if this process is analogous to the 
heating of a compressible fluid flowing in a duct, con- 
sider the case where the duct area is varied as heat is 
added. Let the area vary such that velocity remains 
constant. From the continuity equation, it is seen that 
density must decrease if the product p AV is to remain 
constant. If velocity is constant, the energy equation 
states that the quantity (T q “T) must be constant. Since 
beat is being added and T Q is Increasing, T must be in- 
creasing at the same rate. 

Now consider the incompressible situation. If the 
analogue is to hold, the parameters analogous T and p 
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must increase and decrease, respectively. This is impossible 
since both are analogued by y. 

4. Description of the Water Table. 

The water table used was initially constructed by Naval 
Ordnance Test Station, China Lake, California The design was 
based upon the experience of the General Electric Company. [5] 

Constructed of wood and painted with an acid resistant 
laboratory paint, the table measures approximately 42 inches 
wide by 69 inches long. The legs are fitted with jack screws 
to facilitate table slope changes. 

The original table bed was replaced by a sheet of lucite 
one inch thick, 36 inches wide and 48 inches long. The under- 
side of the lucite was frosted and a two inch grid drawn upon 
it. Three four foot flourscent light tubes were mounted under 
the table in a direction parallel to the flow. The resulting 
illuminated grid is suitable for photography. The sides of 
the channel are made of five-eighth inch lucite and are bonded 
to the table to insure a water tight fit. An adjustable sluice 
was constructed in order to regulate y Q at any given flow rate. 
The bottom edge of the sluice was milled to a one-half inch 
radius to prevent tripping to turbulent flow. 

Water flows under the sluice and across the table from a 
reservoir of three cubic feet capacity. The reservoir is fed 
by a one inch pipe perforated at random intervals to insure 
an even distribution of incoming water. Two screens were fab- 
ricated from honeycomb stock (one-eighth inch cell). These 
were inserted in the reservoir and provide complete "stillness" 
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which insures adequate stagnation conditions. Water enters 
the channel over a lip which has an eight inch radius to pro 
mote gentle transition from zero velocity. 

The water is collected in a return trough and fed back 
to the reservoir through a centrifugal pump and a flowrotor. 
The maximum flow rate of the system is 30 gallons per minute 

Depth measurements are made by means of a three inch 
depth micrometer with a pointed tip attached. The micro- 
meter is bolted to an aluminum bar which extends across the 
table. The bar rests on two steel straps mounted parallel 
to the channel bed. Measurements made with this apparatus 
are accurate to plus or minus 0.005 Inches. This accuracy 
is limited by the bow inherent in the lucity bed and by the 
alignment of the steel straps. Figure 6 shows the table and 
measuring apparatus. 

Detergent was added to the water to reduce surface ten- 
sion and oapillary waves. It was discovered that if enough 
detergent were added the action of the pump caused a suspen- 
sion of bubbles. The resulting flow had a gray color and 
provided excellent contrast for photographs. 

An intravenous feeding apparatus was obtained and used 
to admit drops of dye to the system. Red food coloring (un- 
diluted) proved to be a satisfactory dye. This arrangement 
could be closely regulated , permitting drops to fall on the 
flow at one-half second intervals. Another method of measur 
ing drop interval was to time a number of drops <, say twenty * 
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wiih a stop watch. The drops were then photographed down- 
stream (against the grid) giving an accurate measurement 
of velocity. The fact that the head of the dye reservoir 
is constantly decreasing is insignificant in light of the 
small amount of dyo added ao cue time. This is limited 
by the ultimate coloration of the water . 

5„ 2xanpl.es of Capabilities of the Table 

Photograph s were made of flow about various models 
at several Froude numbers to check the effectiveness of 
the lighting, dye addition system , and constrast. 

Figure 7 shows drops of dye downstream. Twenty drops 
of dye were added over a period cf 12 seconds. 




VELOCITY MEASUREMENT BY DYE DROP METHOD 
18 



The photograph shews the drop 13*2 inches (6.6 grid 
squares) apart. Thus the velocity of the stream is 
1.8 feet per second. Depth measurements taken at the 
same time indicated the critical velocity was 0.8 
feet per second. The Froude number is 2.2. 




TRAILING VORTICES 



Figure 3 ' j a pnotograph cf the vortices prodi d 
K s'ymmetri al airfoil at a sligr.t angle of attack. 
Tn Froud ■ ; . *• ^.6c* Sc } >.as ailovjf-d K i r 

r.rough the pamp tc produ -e visible vortices. 
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FIGURE 9 

STREAMLINE PATTERN PAST A SYMMETRICAL AIRFOIL 



The streamline pattern about the same airfoil is 
sho>m in Figure 9 • Red food coloring was used as a dye. 

The flow past a supersonic airfoil is compared to 
that past a lowdrag subsonic airfoil (of approximately 
equal thickness ratio) in Figure 10. Rote the detached 
"shock" on the plurt airfoil compared to the attached 
shock on the wedge shaped airfoil. The Froude number 
wa 3 ? . 3 • 

Figure 10 shows that oblique shock formed by the 
eenterbody of a ramjet intake at F = 3 . 1 . Note that 
design cercitions are being met; the shock is on the 
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lip of the inleto secondary shocKs are seen in the 
diffuser. 




FIGURE 10 

"SUPERSONIC" FLOW PAST TWO DIFFERENT AIRFOILS 



The hydraulic analog} has also seen used to study 
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cylinder and a converging-diverging nozzle [3] 9 as well as 

tne flow about various ether aerodynamic shapes [10] • Re= 
ferences [7 ] and [10] are considered extremely useful in de = 
monstrating the applicability of the analogy ana the techniques 
required in its use. 

5. Recommendations 

On the basis of knowledge gained during this investigate 
ing, it^is recommended that a new water table be constructed 
for use at the United States Naval Postgraauate School. The 
new table should incorporate a glass bed since the lucite in 
the present table is easily scratched and not easily polished. 
It should also be much larger 9 with the dimensions of the table 
described in reference [i] considered minimum. Other features 
considered mandatory are a variable porosity tail gate and a 
larger settling chamber. The table should also provide for 
application of the “hydrogen bubble “ technique [6] . This 
would eliminate the necessity for a dye addition system;, as 
the bubbles will give accurate velocity and streamline informa= 
tion. 
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